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Abstract: The character of energetic particle transport in the distant heliosheath and especially in the vicinity of the 
heliopause could be quite distinct from the other regions of the heliosphere. The magnetic field structure is dominated by a 
tightly wrapped oscillating heliospheric current sheet which is transported to higher latitudes by the nonradial heliosheath 
flows. Both Voyagers have, or are expected to enter a region dominated by the sectored field formed during the preceding 
solar maximum. As the plasma flow slows down on approach to the heliopause, the distance between the folds of the 
current sheet decreases to the point where it becomes comparable to the cyclotron radius of an energetic ion, such as 
a galactic cosmic ray. Then, a charged particle can effectively drift across a stack of magnetic sectors with a speed 
comparable with the particle’s velocity. Cosmic rays should also be able to efficiently diffuse across the mean magnetic 
field if the distance between sector boundaries varies. The region of the heliopause could thus be much more permeable to 
cosmic rays than was previously thought. This new transport proposed mechanism could explain the very high intensities 
(approaching the model interstellar values) of galactic cosmic rays measured by Voyager 1 during 2010-2011. 
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1 Introduction 

The years 2009-2011 marked the solar minimum in the 
outer heliosphere. During this period intensities of galactic 
cosmic-ray protons and alpha particles measured by Voy- 
ager 1 rose to record high levels. At the time of this writ- 
ing Voyager 1 is at 117 AU heliocentric distance, which 
is at least 23 AU beyond the termination shock. There is 
evidence from plasma flow velocity decrease [5] that the 
spacecraft may be close to the heliopause and is observing 
a nearly pristine interstellar flux of galactic ions. This in- 
terpretation, however, is not consistent with models of the 
global heliosphere that predict a heliosheath thickness of 
~45 AU in the nose direction [9], 

Galactic cosmic ray intensity was expected to be attenuated 
by a wall of enhanced magnetic field near the heliopause 
[8, 7], Magnetic field magnitudes measured by Voyager 
1 up to 2010 [2] are not consistent with the values pre- 
dicted for the magnetic wall [1], which could indicate that 
the spacecraft is still far away from the heliopause. Below 
we show that the magnetic wall is actually quite permeable 
to cosmic rays in the region dominated by the tightly folded 
neutral sheet. The new transport process, first discussed in 
[6] provides a natural solution to the problem of cosmic ray 
intensity enhancement in the distant heliosheath. 


2 Magnetic field structure in the heliosheath 

Because the radial flow velocity vanishes on approach to 
the heliopause, plasma has a long (several decades) resi- 
dence time in the heliosheath. The heliosheath is charac- 
terized by nonradial flows as the plasma is deflected away 
from the radial direction by the obstacle (the heliopause). 
The heliospheric current sheet is stretched in the latitudi- 
nal direction. The process is illustrated in Figure 1, which 
shows the shape of the current sheet projected onto the 
meridional plane. This result was obtained using tracers 
embedded in a 3D flow structure produced from an MHD- 
neutral model of the global heliosphere [6]. The initial 
tilt angle at the Sun was 20° in this simulation. One can 
clearly see an increase in the latitudinal extent of the cur- 
rent sheet in the heliosheath and a compression of the mag- 
netic sectors. Similar results have been previously obtained 
by [4, 1], 

In 2010, Voyager 1 was above the northward extent of the 
current sheet in the supersonic solar wind, but most likely 
well inside the extended sheet structure shown in Figure 
1 [2], It was shown [6] that energetic charged particles 
acquire a degree of mobility across the stack of magnetic 
sectors once the width of the sector d becomes less than 
2r g , where r g is the gyroradius. For a 160 MeV proton, 
this condition means a sector width of less than 0.25 AU 
assuming a background field intensity of |B| = 0.1 nT. In 
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Figure 1: A projection of the heliospheric current sheet 
onto the meridional plane. The termination shock is shown 
with a dashed line (from [6]). 


the solar wind sectors are about 3 AU apart and their width 
decreases to under 1 AU downstream of the termination 
shock. The width continues to decrease into the heliosheath 
without a comparable increase in the magnitude of B (the 
excess plasma escapes between the current sheet to higher 
latitudes). It is likely that the condition d < 2 r g is sat- 
isfied at the current position of Voyager 1. As we show in 
the following section, cosmic ray ions experience enhanced 
transport across a stack of closely spaced sectors. This is an 
essentially kinetic process that depends on a particle’s gy- 
rophase. It is not captured with a nearly isotropic (Parker) 
model of cosmic ray transport [10], 


3 Cosmic ray drift and diffusion in a sec- 
tored field 

When the condition d < 2 r g is satisfied, a charged particle 
is unable to close its gyroorbit inside one sector and will 
cross into an adjacent sector of opposite magnetic polarity. 
The guiding center then undergoes drift motion both paral- 
lel to the current sheets and across the sectors. When the 
sectors all have the same width d, the drift velocity in the 
perpendicular direction is given by [6] 

vdeB 

Vdx / , , \ 5 ( 1 ) 

PC(7T -01-02) 

where e is the particle’s charge, v velocity, p momentum, 
and 0i, 02 are the gyrophases corresponding to the tra- 
jectory crossing points with two adjacent current sheets. 
The drift is bi-directional, i.e., particles can cross the stack 
of sectors both left-to-right and right-to-left, depending on 
their initial position and gyrophase. 



Figure 2: Cross-sector diffusion coefficient k xx as a func- 
tion of sector width relative to the gyroradius. The normal- 
ization chosen is r g = 1 and v = 1. Dotted lines corre- 
spond to k xx ~ D~ x . 


Magnetic sectors observed by the Voyagers in the solar 
wind have variable widths with about a factor of 10 dif- 
ference between the shortest and the longest sectors [3], 
When the sector width varies in a random way, cosmic rays 
undergo random walk (diffusion) in the cross-sector direc- 
tion. Particles reverse the direction of drift motion when 
the next sector encountered is wider than their gyroradius. 
Suppose the sectors are, on average, parallel to the yz plane 
with the magnetic field B = ±By. The trajectory of the 
particle’s guiding center £(i ) in the direction perpendicular 
to the sectors (x) is described by the following relation: 
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where the index n refers to nth sector. The turn-around 
condition is d„ — > r g . 


Guiding center motion is strictly diffusive if d n is a random 
variable. In [6] sector widths were uniformly distributed 
between 0 and D. A trajectory calculated according to (2)- 
(3) has no memory of sectors encountered previously (be- 
fore a trajectory reversal). A particle’s guiding center drifts 
along the sector boundaries (the z direction), which implies 
that the sector width must also vary randomly in the z di- 
rection (i.e., with latitude). 


Here we extend the analysis of [6] with a series of guiding 
center trajectory simulations to reveal the dependence of 
the current-sheet diffusion coefficient on sector width and 
particle’s gyroradius. Figure 2 plots the dimensionless dif- 
fusion coefficient k xx , calculated by fitting the mean square 
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Figure 3: Cross-sector diffusion coefficient k xx as a func- 
tion particle’s speed. Here vo corresponds to r g = D. Dot- 
ted line is k xx ~ v 3 . 

displacement ((x — {x)) 2 ) with a power law with index 1 in 
time, as a function of d/r g (for a fixed r g ). The dependence 
is asymptotically n xx ~ D _1 at large and small D/r g with 
a rapid transition near D/r g = 1. Evidently, when the 
sector width is large, particles tend to drift along a single 
current sheet and their cross-sector transport is suppressed. 

The second set of simulations reveals the dependence of 
diffusion on particle’s velocity. The result for k xx (v)\s 
shown in Figure 3. Here Vo corresponds to a gyroradius 
r g = D. Velocity dependence is well approximated with a 
cubic fit at large and small v with a more rapid transition 
in between. Using B = 0.1 nT and D = 0.2 AU yields a 
perpendicular diffusion coefficient k xx = 1.5 x 10 22 cm 2 /s 
for a 160 MeV proton. This value is well in excess of the 
perpendicular diffusion coefficient predicted by the NLGC 
theory [13]. We conclude that cross-field diffusion in the 
tightly compressed sector field in the heliosheath is very 
efficient which could explain the evident lack of cosmic- 
ray attenuation in the distant heliosheath. 

4 The effect of solar cycle in the distant he- 
liosheath 

Owing to long plasma residence times in the heliosheath, 
a record of the entire previous solar cycle is preserved in 
a compressed form [11, 12], Using the tracer method, we 
performed a simulation of the current sheet with a variable 
tilt angle a given by following expression 

f 2nt\ 

a = 7 + 36.5 ( 1 + sin — J , (4) 

where T = 11 yrs and the angles are in degrees. Figure 
4 shows the minimum (southward) and maximum (north- 
ward) extent of the heliospheric current sheet during a time 
near a negative solar minimum in the inner heliosphere. 


heliographic coordinates 



Figure 4: Regions of different magnetic polarities at the 
time of a negative solar minimum in the inner heliosphere. 
The projections of the termination shock and heliopause 
on the meridional plane are shown with dashed lines. The 
dotted line is Voyager 1 trajectory. 

The two lobes in the north and south are regions of uni- 
form polarity (toward the Sun in the north and away in the 
south). The near-equatorial region in between is filled with 
a warped current sheet and has mixed magnetic polarity. 
This regions extends to high latitudes in the heliosheath (a) 
because of the increase in the tilt angle toward the previous 
solar maximum and (b) because the warped current sheet 
is carried to higher latitudes by latitudinal flows in the he- 
liosheath. 

The current sheet tilt was low during the previous positive 
solar minimum around 1996. Solar wind plasma ejected 
by the Sun during that minimum is still traversing the he- 
liosheath. This region is represented by two thin slices de- 
scending toward the ecliptic from high latitudes in the north 
and in the south in Figure 4. Voyager 1 trajectory passes 
above the lowest latitude of the 1996 solar minimum re- 
gion. Given that the plasma velocity is near zero at Voy- 
ager 1 [5], the spacecraft may well be sampling parcels of 
solar wind emitted at increasingly earlier times. It is even 
conceivable that Voyager 1 will revisit the previous solar 
minimum in the heliosheath before its encounter with the 
heliopause. 

5 Discussion 

Conditions in the heliosheath vary on slower time scales 
than in the supersonic solar wind. One cannot easily char- 
acterize these conditions as specific to the solar minimum 
or solar maximum because both are present at the same 
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time. Voyager 1 may now be sampling solar wind emitted 
by the Sun during the descending phase of solar cycle 23. 
In the framework of this model, one can even interpret the 
long positive sector observed by Voyager 1 at the end of 
2009 [2] as an encounter with the positive unipolar region 
shown in Figure 4. 

The tightly folded heliospheric current sheet is expected to 
dominate cosmic-ray transport near the heliopause, at least 
at low and mid latitudes. As the magnetic sectors become 
thinner, particles acquire mobility in the perpendicular di- 
rection allowing them to cross the modulation wall with rel- 
ative ease. This process could be responsible for the record 
high galactic cosmic ray intensities measured by the Voy- 
ager probes during the past two years. 
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operative agreement with NASA Marshall Space Flight 
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